Width of phonon states on defects of various dimensions 



L. A. Falkovsky * 
Landau Institute for Theoretical Physics, 117337 Moscow, Russia 



Abstract 

We consider the localized phonon states created by defects of various ge- 
ometries near the edge of an optical-phonon branch. The averaged Green's 
function is calculated to study the Raman line shape. The phonon scattering 
by the defects induces broadening and line shape asymmetry. The contribu- 
tion of localized states to Raman spectra has a form of shoulder with a width 
proportional to the square root of defect concentration. 
PACS: 63.20.-e, 78.30.-j 

I. INTRODUCTION 

The interaction of elementary excitations with defects is an old problem of solid state 
physics. For the case of optical phonons, it has attracted considerable attention in recent 
years because of the Raman scattering applications for determination of the isotopic com- 
position ||T|, for detection of stacking faults in polytypes [0,0], for investigations of residual 
strain and strain relaxation at interfaces in semiconductors j|. 

A refinement of these effects takes the interactions among phonons and with defects into 
account. 

(i) The phonon decay into two (or more) acoustic or optical phonons is possible even 
at zero temperature. These unharmonic processes give rise to the so-called natural width 



r nat ~ Uoym/M, where ujq is the optical-phonon frequency, m and M are the electron and 



atom masses. In spite of the small value Jm/M ~ 10 2 , the decay processes give the observ- 



able linewidth in Raman spectra and the quite small mean- free path r 7 ~ ayLu /T nat ~ 10a, 
where a is the lattice constant. The decay cross section depends weakly on the phonon fre- 
quency since no singularities exist in the final density of states for that processes. Therefore, 
the natural optical-phonon width r nat can be considered for given material as constant up 
to room temperatures and for the phonon frequency varying in a small interval of the order 
of r nat . Then we obtain for perfect crystals the Lorentzian line shape of Raman spectra. 

(ii) The phonon interaction with defects should be considered in two aspects. First, it 
gives rise to elastic phonon scattering and makes an additional contribution to the phonon 
width T. This contribution (linear in defect concentration c) can be evaluated || using 
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the golden rule of quantum mechanics. The scattering by defects, being elastic, conserves 
singularities of the phonon density of states. This is important for the Raman scattering, 
because we are interested in phonons with small wave vectors. For the short-range potential 
of defects with the radius of the order of a < r 7 , the singularity arises also in the width 
T as function of phonon frequency u at the phonon branch edge, i.e., at the maximum (or 
minimum) value ujq. Thus we obtain the asymmetrical Raman line shape, that was observed 
in experiments with 3C-SiC H. This asymmetry is sensitive to the defect geometry, which 
may have the point, line or plain form ||. Secondly, the interaction with defects can result 
in appearance of new modes. For instance, the mode forbidden by selection rules in perfect 
crystals at given experimental conditions has been observed in crystals with the stacking 
faults which give an example of plain defects in SiC polytypes 0. Some additional modes 
were found also in computer simulations for the sequence of number layers with the stacking 
faults ||. In addition, the localized mode can be associated with defects 0. 

The problem of localized states on the single isotopic defect has the exact solution. 
However, the case of the low, but finite defect concentration involves serious difficulties 
||||. That is because the solution should be self-consistent for the large distances essential 
near the edge of phonon branches. If the natural width is ignored, the answer obtained in 
the series expansion of the phonon self-energy is the following. The edge of phonon branches 
is displaced linearly in c and the band of localized states arises with a width of the order of 
c 2 . This c 2 behavior is explained by fluctuations in the spatial distribution of defects. 

The aim of the present paper is to consider (taking the natural width into account) the 
situation when the states localized on defects can exist. We are especially interested in the 
case when the localized states appear near an edge of optical-phonon branches. Broadening 
of localized states is found just in the lowest order in the defect concentration c, if this addi- 
tional width is larger than the natural phonon width. This unusual behavior is conditioned 
by the resonance interaction of localized states with the close edge of the extended phonon 
branch. The defects of different geometries are treated such as point defects, line defects 
involved by dislocations, and plain defects produced by stacking faults or by the crystallite 
boundaries. They give the different broadening and asymmetry of the Raman line shape. 
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II. INTERACTION OF PHONON MODES WITH DEFECTS 

We solve the equation for the phonon Green's function in the form: 

D(r, r', u) =D {r- r', u) - J dv"D (v - r", u)U(r")D(r", r', cu), (1) 

where the interaction with defects located in points r n is 

U{r) = $>(r - r n ) = £ W (q)e^ r - r «). 

n n,q 

For an isotopic defect, u(q) is constant: w(q) = u = (M — M)a 3 /M , where M and M 
are the masses of the defect atom and the host atom, respectively, a is the lattice constant. In 
the following, the large distances from the defect on the atomic scale are important. Then 
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for any type of defects, we can set the Fourier component w(q) equal to the constant uq 
supposing the short-range potential. For the same reason we expand the phonon spectrum 
near the Brillouin-zone center taking the Green's function in the absence of defects as 

D (k, uj) = {uj 2 - s 2 k 2 -uj 2 - iujr nat )-\ 

where the natural phonon width r nat is connected to decay into acoustic or optical phonons. 

The imaginary part of -D(k, uj) gives the differential density of squared frequencies. In 
addition, the Raman intensity I{oj) for phonon excitation with frequency uj and momentum 
k is proportional to InxD(k, uj) which gives the Lorentzian line shape in the absence of defects 

^pnat 

I(uj) ~ lHLD (k, uj) 



uj 2 q - s 2 k 2 - to 2 ) 2 + (wr nat ) 2 ' 



if T nat is assumed to be small in comparison to the frequency uj . 

In the presence of defects, the Green's function must be averaged over defect positions. 
It is convenient to use the well-known diagram technique for impurities scattering |TT| . This 



method can be modified [12] in order to take the possible localized states into account. 



Performing the iterations of Eq. ([!]) in U(r) and averaging over the defect positions, we 
have for the sum of noncrossing diagrams: 

D(k,u) = D (k,u) +cD 2 (k,u) (2) 

00/ \ n 

xE(-«o)" +1 (cD (k,oj) + Y,D(q,u;)) , 

n=0 V q I 

which gives after some algebra 



D(k, uj)' 1 = D (k, uj)- 1 + cu (l + u E D(q, u)^ 



(3) 



The term of first order only in the defect concentration c is held in the phonon self-energy 
given by the second term on the right-hand side of Eq. ([3]), since we assume that the distance 
between defects (r c oc c" 1 ^ 2 for line defects) is larger, than the phonon mean-free path r 7 . 

Let us emphasize that the phonon self-energy involves the perturbed function -D(q, a;), 
but not .Do(q, a;). Then the Green's function -D(q, uj) has to be found from the integral 
equation Integration (summation) is perfomed with respect to the 3D- vector q for 

point defects, over the 2D- vector q^ for line defects (then q z = k z ), and over q z for plain 
defects (then qy = ky). 

The integral equation (|3|) can be reduced to an algebraic equation. We introduce the 
new unknown function ( and the complex variable Co in the case of point defects by the 
definition 

D(k, u) = (C - s 2 k 2 Y\ ( = uj 2 -uj 2 - iujT nat . 

Then Eq. (|J) gives 



Co + co; 3 (^-«: + ^-C) • (I) 
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In the case of line defects 



D(k,u) = (C-s 2 kl)- 1 , 

and we obtain 

C = Co + cul 

For plain defects 

D(k,co) = (C-sXr 1 , 

we have 

C = Co + cul 

where A is the dimensionless coupling constant proportional to uq, c is the atomic concen- 
tration of defects in 3D, 2D, or ID region, respectively for the point, line, or plain defects, 
k is the cutoff frequency, the z axis is taken in direction of the line defect or normal to the 
plane defect, kii is parallel to the plane defect. The branches of yj — C and In (— C) have to 
be taken in the upper half-plane. We consider the retarded Green's function which has no 
singularities in the upper half-plane of complex variable u>. 

The equations have two singularities. One at C = is related to the branch 

edge of the extended optical-phonon mode. Other at (i is related to the localized mode and 
determined by the zero of expressions in parentheses of Eqs. (f|)-(0). For instance, for line 
defects we have 



which gives 

= -K 2 e' 1/X . 

Notice, that the zero Cz of the expressions in parentheses, if it exists, is negative. The zero 
can exist, if A > 1 in the case of point defects, and for any positive A in the cases of both 
line and plain defects. 

The situation is schematically explained in Fig. 1, where Co determined by Eqs. (f|)-(§P 
is shown as function of C- In the regions plotted by solid lines and indicated by the arrows 
1, the solution ( is real for real values of Co and has no singularities in the upper half-plain 
of u, as it must be for the retarded Green's function. If r nat = 0, this region makes no 
contribution to the density of phonon states. The second solution ( (dashed lines indicated 
by the arrows 2) corresponding to the same Co should be rejected since it has singularities in 
the upper half-plain. For real Co ; the solution £ has the imaginary part in the regions shown 
by dot-dashed lines and indicated by arrows 3 and 4. The region 3 represents the extended 
phonon mode. The region 4 shows the localized states. 

a) Low concentration of defects: unrenormalized approximation. 

We can solve Eqs. analytically supposing c to be small. Using the iterative 

method (if |C — Col "C |Co|) we obtain on the first step for point defects 



2 2/2 2 • -piiat 

Co = UJ Q — S k z — UJ — iujY 



(5) 



/- 2 2 ; 2 2 • -nr 

= UJ Q — S — UJ — tUJL 



(6) 
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-1 



C = Co + CW - - K 



(7) 



for line defects 



c 



Co + cul 



111 



-Co 
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) 



(8) 



and for plain defects 




) 



-i 



c 



%/=Co~ 



(9) 



This unrenormalized approach corresponds to substituting -D (q, u) for D(q,u) in the self- 
energy of Eq. (^) . 

In Eqs. (0) - ©) the terms proportional to the concentration c represent the contribution 
of defects. If Re Co > 0, i.e., for the region shown by the arrow 3 in Fig. 1, the frequency 
belongs to the extended states. Then the imaginary and real parts of the second terms on 
the right-hand side of Eqs. (0)-© give the phonon shift and the additional width induced 
by the defect - phonon scattering. If Re Co < and for positive A, the localized states can 
exist. In the absence of the natural width, they always exist for the cases both the line 
and plain defects, in accord with the well-known statement of quantum mechanics about 
the bound state in a shallow potential. However, for the finite value of r nat , the coupling 
constant should be strong to give the substantial gap: W| - Uo > Y na,t . This condition can 
be rewritten as ri <C r 7 , if we introduce the radius of localized states ri = aJtu / (uoi — loq) 



defects c. It is hard to satisfy the conditions of validity, if Co — ► Ci- I n fact, the equation 
Co = (i determines the frequency ui of localized phonon states. For |Co — Gl ^ |Cz|> we 
expand the expressions in parentheses of Eqs. (§])-© m (C — Ci) an d arrive to the quadratic 
algebraic equation which has the solution 



where b = 4(— QY^Uq/t: for point defects, b = —Q^o for line defects, and b = 2tu (—Q) 3 ^ 2 
for plain defects. We notice again that the solution with the positive imaginary part should 
be taken in Eq. (|10|). 

If Co — > Ch we obtain a negative value under the square root in Eq. (|r0|), that gives the 
imaginary part of C- Then the imaginary part is determined in competition of u Q r nat with 
2 (6c) 1 / 2 . The region where T nat < 2{bc) 1 ^ 2 / ujq is indicated by the arrow 4 in Fig. 1. This 
condition may be compatible with the conditions (i) r c ^> r 7 for the expansion in terms of 
c used in Eq. ([|) and (ii) r 7 3> r\ for the existence of localized states. Here, the localized 
states make the maximal contribution to the differential density of states and to the Raman 
intensity, 





(10) 
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which is proportional to the square root of the defect concentration. In the opposite limit 
pnat > 2(6c) 1 / 2 /co , o, we can use the Eqs. (0)-© which give the contribution linear in the 
defect concentration. 

c) Self- consistent approach: the numerical solution. 

As a matter of fact, each of the equations (||), (H), and (Q) is a system of two equations for 
the real and imaginary parts of (. This system may be solved numerically. It is convenient 
to introduce the unknown real variables x and y as 



supposing y > and —ir/2 < arctan z < 7r/2 in order to comply with the analytical 
requirements in the upper-half plain. The similar form is used also for £ with x — > uo 2 — uj^ 
and y — > c<jr nat . 

The results of the numerical solution are shown by solid lines in Fig. 2 (line defects) and in 
Fig. 3 (plain defects) for the cases, when the coupling constant A is negative and no localized 
states exist. We used the values of c^o — 520 cm -1 and r nat = 3.2 cm -1 , which correspond to 
pure silicon, and the cutoff parameter k = ujq. The unrenormalized approximation is shown 
by dashed lines. We see that the interaction of phonons with defects results in broadening 
and shift of the Raman line from the position (dotted line) in the absence of defects. 

The broadening is much more on the low-frequency side of the line than on the high- 
frequency side, that gives the line shape asymmetry. This is an effect of the phonon density 
of states: the elastic interaction with defects makes a contribution to the phonon life-time 
only in the region uo 2 < u 2 ., where the final phonon states exist. 

In Fig. 4, the dependences of the line shift and width on the concentration are shown 
for the point, line, and plain defects. The value of coupling constant A = — 1 is taken. We 
define the full width at half maximum F (FWHM) in the usual way: T = uj + — cj_, where 
the Raman intensities I(u±) = 0.5/(u; max ) and u + > a; max > oj_. The behavior of the width 
T and of the line position cu max , calculated with the help of Eqs. (f|), (|j), and @ is shown in 
the left panel. In the right panel, the Raman line asymmetry {oj + — co> max ) — (c<-> max — ujJ) is 
plotted. We see that the plain defect affect more intensively on the width and asymmetry, 
than the line and point defects. To the contrary, the Raman line is shifted more for the 
point defects. 

The cases of the positive coupling constant, when the localized states arise, are shown 
in Figs. 5 and 6. The accurate results (solid lines) differ considerably from that obtained 
in the unrenormalized approximation (dashed lines). The maximum that was found in the 
unrenormalized approximation is smoothed out, transforming into a shoulder with the width 
of the order of a/c in accord with Eq. fllQl). Emphasize, that the total number of localized 
states, integrated over frequency, is proportional to c independently of the approximation. 

The difference in results for defects of different geometries is evident from Figs. 2- 
6. These distinctions in the Raman line shape and shift can be detected in fitting to the 
experimental data. 
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III. CONCLUSIONS 



To summarize, we investigate here the structure of the phonon density of states and 
the Raman line shape in the situation when localized levels u)\ interact strongly with the 
closely-spaced edge of the extended optical-phonon states. In the linear approximation in 
defect concentration c, the contribution of localized states does not have the well-known 
form of ^-function: c6(u — u>i), if the effect of defect concentration is more essential than the 
phonon natural width. 

The line shape of bulk-phonon extended states is asymmetrical. The additional width 
(proportional to the defect concentration) arises on the side of line related to the band of 
extended states. As a most intriguing example, we would like to compare the behavior of the 
Raman lines of graphite at 1600 cm -1 and of diamond-like material at 1300 cm" 1 . When 
defects are added to graphite, the low-frequency wing of the line becomes to drop more 
slowly than the high-frequency wing (see Figs. 1 and 2 in Ref. [O]) in just the same way as 



shown in Figs. 2 and 3 of the present paper. But if it is diamond samples that are impaired, 
the opposite, i.e. high-frequency, side slopes more gently (see Fig. 7 in Ref. Jl3| and Fig. 2 in 
Ref. [0). Using the present theory, the explanation of the phenomenon may be attributed 
to the difference in the phonon dispersion of these substances. Namely, contrary to graphite, 
there is a minimum of the optical-phonon branch at the Brillouin zone-center in diamond 
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FIG. 1. Plot of function Co(C) given by Eqs. dH|). In the regions shown by the dot-dashed 
lines (indicated by arrows 3 and 4), the solution C(Co) has a finite imaginary parts for real Co- The 
region 3 corresponds to the extended phonon mode. The region 4 is related to the localized states. 



9 



4.5 



3.5 



CO 

■E 3 



^2.5 

CO 
CD 

CO 



CO 



1.5 



0.5 







I 


I I I 

. ;^X0.75 


- 


c = 


- 


line defects 


1 

1 

1 

1 

1 


X = -1 

[ : 
l : 
l ; 


1 

1 

- c = 0.1->/ 


I : 
I : 

:'l ; . 

: I : 
: I : 


■ H- ■ ' ' 


V . 

i ^ — i 



480 



580 



500 520 540 560 
frequency transfer (1/cm) 

FIG. 2. Imaginary part of the phonon Green's function versus frequency to for parameters of 
silicon: loq = 520 cm -1 , r nat = 3.2 cm -1 . The number of line defect per atomic volume c = 0.1 
and the phonon-defect coupling constant A = — 1 (the solid line for the self-consistent solution, the 
dashed line for unrenormalized approximation). The results for pure crystal (c = 0) are shown by 
the dotted line. 



10 



4.5 



3.5 



C/i 

E 3 



CO 



2.5 



CO 
CD 



CO 

E 

CO 

DC 



1.5 - 



0.5 







I 


. : ^x0.75 


_ 


c = 


- 






i plain defects 


I 


( r 

t r A, = -0.4 

1; |i 

: y 

: ( : 


I 

I 

- c = 0.1 ->/ . 


1: 

l' 
1 

i: 


i ■ 


i — — i i 



480 



580 



500 520 540 560 
frequency transfer (1/cm) 

FIG. 3. Same as Fig. 2 but for plain defects with the coupling constant A = 



-0.4. 
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FIG. 4. The theoretical Raman line width (upper part of left panel), shift (bottom part of left 
panel), and asymmetry (right panel) as functions of the defect concentration. The phonon-defect 
coupling constant A = — 1. Results for the point, line, and plain defects are plotted by the solid, 
dashed, and dot-dashed lines, respectively. 
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FIG. 5. Same as Fig. 2 but for the positive coupling constant A = 0.4, when localized states 
exist. The defect concentration c = 0.02. 
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FIG. 6. Same as Fig. 2 but for plain defects with the positive coupling constant A = 0.4 and 
the concentration c = 0.04. 



14 



